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The decays of r ^ Sttz/ and r -^ nK*^, Kpu are calculated using the hard 
^ i pion and kaon current algebra and assuming the Axial- Vector meson dominance 

^ ! of the hadronic axial currents. Using the experimental data on their masses and 

widths, the r decay branching ratios into these channels are calculated and found 
to be in a reasonable agreement with the experimental data. In particular, using 






Oh! 

^ • the available Aleph data on the 3tv spectrum, we determine the Ai parameters, 

niA = 1.24 ± 0.02GeV, Fa = 0.43 ± 0.02 GeV; the hard current algebra calculation 



X 



yields a 37r branching ratio of 19 ± 3%. 
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Application of hard current algebra in r decay was initiated soon after its 
discovery [1, 2, 3]. It was also pointed out that in a related calculation, using the 
hard current algebra technique, the cross sections for e'^e~ -^ An can be calculated 
and agree with the data in the 1 GeV region within a factor of 4 instead of a factor 
of 10 using the usual soft current algebra. This technique was extended to the ^64 
decays with an implementation of the unitarity in the 2 pion channels [4] , and was 
later used with success in the resolution of the 77 — *> Stt problem [5] . 

The technique of the hard pion current algebra consists in using the PCAC 
and the Lehmann Symanzik Zimmerman reduction formula [6] by contracting out 
the pseudoscalar fields. One would get the expressions for single, double or triple 
equal time commutator relations (ETCR) whose Fourier transforms are the physical 
matrix elements involving soft pion emissions. The remaining terms, which cannot 
be calculated by current algebra technique, involve higher order in the pseudoscalar 
meson momenta and therefore do not contribute in the hypothetical world where 
the pseudoscalars have zero four momenta. In the physical world where their four 
momenta do not vanish, corrections must be taken into account in an approximate 
manner by using substracted dispersion relations; the substraction constants are 
given by the current algebra low energy theorems and the substraction points are 
made at the scales where these current algebra low energy theorems are valid [7]. 

To illustrate our point let us consider the e+e~ -^ An. By contracting out the 
2 S-wave pions we obtain the single and the double ETCR. The single ETCR term 
is related to the physical matrix element r -^ 3nu where the axial vector meson 
Ai dominates; the double ETCR is related to the physical pion form factor which 
is dominated by the vector meson p. This technique enables us to incorporate the 
efi^ect of the heavy fields p and Ai which are the main features of the low energy 
QCD and in a way which is consistent with chiral symmetry and experimental data. 
This method is however not unique. 

Instead of this approach, it is now a fashion to use the effective Chiral La- 
grangian which treats elegantly the low energy theorems involving the Nambu- 
Goldstone bosons; quantum corrections are treated perturbatively; this method is 
effectively a power series expansion in momenta and hence cannot take into account 
of the resonance effect. One is forced to accept the fact that Chiral Perturbation 
Theory, as is usually practiced, described only the low energy tail of the resonance 
which is certainly not the gross feature of the physics involved. The remedy for 
this approach is to build in the theory the heavy fields p and Ai and possibly also 
a heavy a fields as was done by a number of authors [8, 9]. It is suggested here 

1 



that the loop corrections should be treated unsystematically in the bubble chain 
approximation in order to make these heavy fields unstable in a way which is con- 
sistent with the unitarity requirement. This procedure is the same as the usual way 
of handling the W and Z propagators in the standard model calculation. One then 
can incorporate many nice features of the old Vector Meson Dominance models but 
taking into account also of the low energy chiral properties of the pseudoscalars. 
This approach will be explored in the future. 

The purpose of this letter is to study the simpler processes r -^ 3tvu and 
r -^ 7riir*z/, Kpu and leave the more complicated process e~^e~ -^ An^ or r ^ Atiu 
for a future publication. We do not expect to achieve here the precision of the 
order of 10% or better which is usually obtained for soft pion emission processes 
like -ftr/2, Kl^ and -ftr/4 or the S wave ttN scattering lengths etc. This is so because 
the matrix elements depend on the scalar product of pion momenta to that of the 
current which is large in the physical region and considerable correction has to be 
made in order to reach the chiral limit. 

Using the hadronic properties (widths and masses) of the Axial-mesons Ai, Qi 
and Q2 and treat them as unstable particles in the usual way, their branching ratios 
and spectra in r decays are calculated and found to be in good agreement with the 
experimental data. 

We begin first by recalling the well-known formula [10] for the ratio 
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where ao((5^) and ai((5^) are respectively the spin and spin 1 part of the hadronic 
spectral functions and Q^ is the Cabbibo angle. The expression multiplying with 
cos^ dc is for the IS.S = hadronic tau decay and that multiplying with sin^ 9c is 
for /S.S = 1 hadronic tau decay. 

We want to study the matrix element for the axial vector hadronic current 
involving three pseudoscalar mesons. For clarity, we first make the simplified ap- 
proximation that the system of three pseudoscalar mesons can be represented by 
a pseudoscalar and a vector meson. This assumption is reasonable because in the 
usual angular momentum decomposition one pair of the pseudoscalars have to be 
in the relative P state and will be shown to be dominated by the vector meson. We 



assume furthermore that the axial hadronic currents are dominated by the Axial 
vector mesons Ai, Qi and Q2, just the same as the vector hadronic currents are 
dominated by the vector mesons p and K* . A more precise study of r ^ 3nu is 
given, treating the vector meson p as a resonant 27r state. 



I) r ^ Sttz/ Decay 

Let us first consider the AS = decay. As mentionned above we approximate 
the Stt state by a rrp state. The most general matrix element can be written as: 

{n-{k)p^{p)\Al-^'mo) = h{Q^)e, + e.k {{k+p)^MQ') + {k-p),h{Q^)) (2) 

where Q^ = (k+p)'^ and e is the polarisation vector of p. /i, /2, and /s are complex 
form factors and are only functions of Q^. Current algebra soft pion theorem, which 
is obtained by taking the limit k^ -^ 0, gives only information on /i but not on 
the other 2 form factors. In an explicit model, it was shown that they contribute 
little to the r -^ wpu. Interested readers are refered to the original article [5]. (We 
assume here that the decay constant of n' is sufficiently small and hence can be 
neglected). Using the standard low energy current algebra theorem and taking the 
limit /c^ ^ we have: 

]im{n-{k)p'{p)\Al-^'{0)\0) = -V2^e,{p) (3) 

where f^^ = 93MeV, and fp is defined by the rate of p ^ e+e". Using the ex- 
perimental data [11] we obtain, fp = O.llSGeV^. This value of fp is equivalent to 
writing approximatively the pion form factor as 

F7r(s) = ml{l + ds/m'^)/ {w?p — s — inipV p{s)). A good fit to the experimental data 
is obtained with 5 — 0.2. In fact the more general form of Eq(3) reads 

lim (7r-(fc)7r+(g07r-(?2)l4-^'(0)|0) = -^F.(.)(qi - q^), (4) 

For convenience we shall first use Eq(3). The 37r matrix element below the pyr 
threshold can be straightforwardly obtained from Eq(4). Using Eq(2) in (3) we 
have: 

/i(mj) = -v^f (5) 



Let us start with the narrow width approximation for the Ai propagator. Using Ai 
dominance for the form factor we have: 

The generahsation of Eq[6] to take into account of the unstable nature of Ai can be 
straightforwardly made. Using the Ai dominance hypothesis for the axial current, 
the general expression for /i((5^) is: 

f m\ — ml, — TT(mi) 
h{Q') = -V2^ i ; \-0- (7) 

where we use the standard prescription for describing an unstable particle, with 
7r(Q^) being the self energy operator of the Ai resonance and is obtained by the 
bubble summation of the ixp intermediate states, similarly to the treatment of the 
W and Z propagators in the standard model. In order to have the usual Breit 
Wigner description of a resonance, we must make a twice substracted dispersion 
relation with i?e[7r(?B^)] = i?e[7r'(m^)] = [12] where tua is the Ai mass: 



(Q^ - m\Y f Im[n{z)]-Im['n:{m\)]-{z - m\)Im[7i'{m\)] 
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where we define the 7r*^p+A^ vertex as g Aptt^{A) .e{p) , A (Q^jTM^, m^) = (Q^ 
{rUp + m7r)^)((5^ — {rUp — TOtt)^), and P stands for the principal part integration. 
The dispersion integral (8-a) can be written as: 



where 
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-^i(Q^); -^2(Q^) and l2{Q^) are defined in the appendix. Using Eq[7] the spectral 
functions ao(<5^) and ai((5^) can be straightforwardly calculated: 
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We want now to generalize Eq(8-b) and Eq(lO) to take into account of the 
unstable nature of the p meson, i.e we want to treat it as a resonant 2n P state. 
Instead of Eq(8-b) we have now: 
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Instead of Eq(lO) we have: 
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(12) 
where Tp{s) is the p width; in terms of the coupling constant (/pttti- we have 

mpTp{s) = -j^s f a/1 — 4mf/sj . From experimental data on the p width, we 
obtain (/pttti- = 6. 

Re['K{Q'^)] is calculated as in eq[8-a] with Im[K{Q'^)] replaced by Im[K{Q'^)] 
given in Eq[ll]. Re[Tf{Q'^)] is calculated numerically. The total branching ratio for 
Stt is obtained by multiplying Eqs(lO) and (12) by a factor of 2. 



Using the narrow width approximation for the p resonance, the Tip spectrum for 
nfiA = 1.2GeV and Fa =0.35, 0.4 , 0.45 GeV is shown in Fig (1). The corresponding 
R ratios are respectively 2, 1.8, 1.6 which are much larger than the experimental 
data. Because of the p narrow width approximation, we cannot calculate the Stt 
spectrum below the np threshold s = {nip + rriT^)'^. Furthermore the up spectrum 
is found to be innacurate near to the np threshold. 

We show below these disagreements with the experimental data are due to the 
zero width approximation for the p resonance. 

In Fig(2), i?e[7r((5^)] in the Tip approximation, and the more correct form 
-Re[7f((5^)] are shown. It is seen that Re[Ti{Q'^)] has a much smoother behavior than 
the approximate one i?e[7r((5^)] near the Tip threshold. Note that at the current 
algebra low energy theorem limit s = rri^, Re[fi{Q'^)] is smaller than Re[7i{Q'^)]. 

In Fig 3, 4 and 5, the Stt spectra are plotted as a function of the Stt invari- 
ant mass where, for clarity, we normalise the theoretical prediction to the number 
of events. The data are taken from the ALEPH group [15] for illustration pur- 
pose. Because the experimental correction for acceptance has not been made, the 
experimental results should be considered as preliminary. In table 1 the values of 
Rstt = p/^g^^j are given as a function of F^ and uia- 

It is seen that Rst^ is an increasing function of m^ and decreasing function of 

If we assume that the acceptance correction to the ALEPH data was negligible, 
our best fit to the Stt spectrum gives the following ranges of m^ and F^: niA = 
1.24 ± OmGeV, Ta = 0.43 ± 0.02 GeV. 

Using these results, the calculated branching ratio for r -^ 3tiu is 19 ±3%. The 
central value corresponds to our best fit which is shown in fig(6). This value agrees 
with the ALEPH branching ratio 19.14 ± 0.48 ± 0.44%. 

The main difference between our results and the recent studies [16,17] of the 
T — s> 37rz/ is that in our approach we can make more predictions due to the appli- 
cation of the low energy current algebra theorem. That is, using the experimental 
mass and width of the axial vector meson Ai, we can predict the r ^ Aiu branching 
ratio. A more ambitious approach to reduce the Ai parameters is to predict the Ai 
width as a function of its mass, similarly to the KSRF relation to the p meson[18]. 
This calculation will be dealt in a forthcoming paper [19]. 
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II) T -^ Kpy and K*7rz/ Decays 

The study of the /S.S = 1 decays is more comphcated due to the presence of 
the two axial vector mesons (5i(1270) and (52(1400). It is an experimental fact that 
Qi is coupled strongly to the pK and weakly to nK* while Q2 couples strongly 
to TT-ftT* and weakly to Kp. Using the experimental data B.R{Qi -^ Kp) = 42%, 
BR{Qi -^ ttK*) = 16%, B.R{Q2 -^ Kp) = 3% and BR{Q2 -^ nK*) = 94% 
together with the total width Tq-^ = 90MeV, Tq^ = 170MeV we obtain (in unit of 
GeV): \9q,pk\ = 2.6±0.20, \9q,pk\ = 0.62±0.30, \9q,^k*\ = 1.18±0.20,\gQ,^K*\ = 
3.66 ±0.10 

Since Qi and Q2 are a linear combination of the states from two different 
octets, one cannot make use of the SU(3) symmetry to fix the signs of the coupling 
constants g [13]. 

Similarly to Eq(3), we can derive the low energy theorems for the matrix ele- 
ment (p°(p)i^~(/c)|Afj~*^(0)|0) by taking the kaon soft, and the low energy theorem 
for the matrix element {n~(k)K*^{p)\A'lj~^^{0)\0) by taking the pion soft. 

In the narrow width approximation for Qi and Q2 resonances, we would expect 
to have the following low energy theorems: 

f 9QlpK , n 9Q2PK _ fp 



+ 



1 p 2 p J t^^ ('\'\\ 

gQiTTK* f 9Q27rK* _ Ik* 



m{ - mj^, mi - mj^, f^+ 

where /q-^, fq^ are similarly defined as fp. 

Eq(13) are strictly not correct because of the unitarity constraints due to the 
complication of the coupled channel problem. A correct treatment of this problem 
will be the subject of a separate publication. Because the experimental data on 
r — ^ K~ p^u and r -^ ti~ K*^^ are very crude, we make the simple approximation 
that Qi couples only to pK and Q2 to K*tx and therefore neglect the mixing of 
these 2 channels. Solving Eq[9] for fq^ and fq^ we have, fq^ = 0.27GeV^ and 
fq, = 0.39GeV^. 

Similarly to /i(s) defined in Eq(2), instead of Eq(7), we now have : 

fKpiQ' 
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where tti{Q^) and 7r2((5^) satisfy a twice substracted dispersion relation with 
Re[TTi{ml)] = Re[7t[{ml)] = and Re[iV2{ml)] = i?e[7r^(m^)] = 0. It should be 
noted that because Kaons are made of up, down and strange quark, the chiral 
(soft) Kaon limit must be taken simultaneously with the soft pion i.e we must use 
the SU{3)l X SU{3)r limit. Because the phase space for Kp is tiny as compared 
with that of K*tt in Qi decays, we must make a correction in the expression for 
7ri(s) to take into account of this special situation. More explicitly we have: 
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(15) 

It is a good approximation to make a 5 function for the K* propagator on the right 
hand side of the equation (11). Similarly we can make the 5 function approximation 
for the p and K* propagators. 



T \ rn2^l o,n 9Q,pK^{Q'':rn%m\) ( X{Q\ml,m],, 
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where rx*(s) is the K* width. A straightforward calculation gives: 

2 Q 

rriK'^K'-is) = ^^27r^ V'^i^^ ''^'k^ '^^) /'^^ • F^om the experimental data on the K* 
width, we obtain Qk^-kK = 4.48. 

The self energy operators tti and 7r2 can be straightforwardly calculated using 
dispersion relations. The corresponding spectral functions are easily calculated. 
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and a similar expression for oq. Numerical calculation gives: 

Si?(r -^ p^K-p) = 0.1% 



(16) 



BR{t -^ n-K''*u) = 0.4% 

These results are in good agreement with the values of TPC/Two-Gamma collabo- 
ration [U]:B (r -^ K*^n-u, neutrals) = 0.51±0.2±0.13 and B {t ^ R-n+TT-u) = 
0.7 ± 0.2. The K~ p^u mode is therefore consistent with zero. 
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These integrals are conveniently expressed in terms of a generating function 
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For convenience we give the analytic continuation of this function to other 
regions. 
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Ta-- 


= 0.35 




mass (in GeV) 


-Rstt 


1.15 


1.02 


1.175 


1.12 


1.2 


1.20 


1.25 


1.40 


Ta = 0.4 


mass (in GeV) 


-Rstt 


1.15 


0.88 


1.175 


0.96 


1.2 


1.04 


1.26 


1.20 


Ta - 0.45 


mass (in GeV) 


-Rstt 


1.15 


0.76 


1.175 


0.84 


1.2 


0.9 


1.25 


1.05 



Table 1 : Results for Rs-^ defined in (1) for different Ai parameters 
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FIGURE CAPTIONS 

Fig.l Prediction for r -^ 3nu spectrum as a function of the Stt invariant mass 
in the narrow width approximation for the p resonance {2tt P state = p). 
dashed/sohd/dot-dashed curves correspond respectively to r^=0.35, 0.4, 0.45 
GeV, the Ai mass is fixed to 1.2 GeV. The experimental data are those of 
Aleph group [15]. 

Fig. 2 Comparaison between Re[TT{Q'^)] and Re[K{Q'^)] (solid/dahed curves) as de- 
fined in (8-a) and (11) (in units of GeV'^). 

Fig. 3 Calculation of the Stt invariant mass spectrum for r -^ 37tu decay using 
Ta = 0.35GeV. Long-dashed/ dot-dashed/short-dashed/solid curves corre- 
spond respectively to niA =1.15, 1.175, 1.2, 1.26 GeV. 

Fig. 4 Idem using r^=0.4 GeV. 

Fig. 5 Idem using r^=0.45 GeV. 

Fig. 6 Our best fit for the r -^ 3ttu spectrum correspending to niA = 1-24 GeV, 
Ta = 0.43 GeV. 
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